Viscosity control of the dynamic self-assembly in ferromagnetic suspensions 
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Recent studies of dynamic self-assembly in ferromagnetic colloids suspended in liquid-air or liquid- 
liquid interfaces revealed a rich variety of dynamic structures ranging from linear snakes to axisym- 
metric asters, which exhibit novel morphology of the magnetic ordering accompanied by large-scale 
hydrodynamic flows. Based on controlled experiments and first principle theory, we argue that the 
transition from snakes to asters is governed by the viscosity of the suspending liquid where less 
viscous liquids favor snakes and more viscous, asters. By obtaining analytic solutions of the time- 
averaged Navier-Stokes equations, we gain insights into the role of mean hydrodynamic flows and 
an overall balance of forces governing the self-assembly. Our results illustrate that the viscosity can 
be used to control the outcome of the dynamic self-assembly in magnetic colloidal suspensions. 

PACS numbers: 87.19.ru, 81.16.Dn, 75.50.Tt 
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Fundamental principles guiding self-assembly in non- 
equilibrium colloidal systems continues to attract enor- 
mous attention in physics and engineering communities 
jll-[l3| . The interest is stimulated by the need for creating 
smart materials capable of self-assembly, adaptation, and 
for the design of tunable structures that can perform use- 
ful tasks at the microscale 
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including targeted cargo 
delivery [l5| , stirring in microfluidic devices |16j , and con- 



la 



trol of optical properties of the media 

Studies of dynamic self-assembly in ferrornagnetic col- 
loids dispersed at liquid-air interfaces [19|, |20[ and en- 
ergized by an alternating {ac) magnetic field revealed 
highly organized, dynamic linear structures - magnetic 
snakes. The snake emerges spontaneously from a random 
dispersion of particles in a certain range of frequencies 
and amplitudes of the ac magnetic field. While for low 
frequencies of the applied magnetic field the snakes are 
immobile, with the increase in frequency they turn into 
self-propelled entities [21|. Surprisingly, fundamentally 
new structures - localized magnetic asters and arrays of 
asters - emerge when the same colloidal suspension is 
confined at the interface between two immiscible liquids 
and is energized by the alternating magnetic field 22 1. 

Both magnetic snakes and asters generate complex 
fiows in the fluid and possess magnetic ordering and dy- 
namic organization highly unfavorable under equilibrium 
conditions. While magnetic snakes are essentially lin- 
ear and comprised of antiferromagnetically ordered seg- 
ments of ferromagnetically ordered chains of micropar- 
ticles [13, asters develop radial structural order with 
the ferromagnetically ordered chains emanating from the 
center of each aster |22|. The mean fiows excited by 



the snakes and asters have fundamentally different mor- 
phology; snakes create quasi-two dimensional fiows with 
quadrupole symmetry confined near the surface 120| and 
asters induce three-dimensional toroidal bulk flows 22 1. 



The main forces that control dynamic self-assembly in 



such systems involve not only magnetic dipole-dipole and 
steric interactions between the particles but also non- 
trivial hydrodynamic forces stemming from deformation 
of the interface, viscous drag, and entrainment by the 
large-scale mean flow. The striking difference between 
self-assembled structures in liquid-air and liquid-liquid 
systems remained unclear since both systems were driven 
similarly. Thus, it is critical to understand the funda- 
mental physical parameters controlling the transition be- 
tween these two distinctive dynamic states. 

In this Letter, we perform a systematic experimental 
and theoretical study of the snake-aster transition. It is 
widely believed that because the motion of each individ- 
ual colloidal particle is strongly ovcrdamped, the viscos- 
ity sets only an overall time scale. Moreover, the motion 
of fluid is often described by the linear Stokes equation, 
as, e.g., in Ref. [23|, and admits a one-way coupling be- 
tween the solvent and the particles, when only the par- 
ticle dynamics is influenced by the liquid flow but not 
vice versa [2J| . For our system, neither of these assump- 
tions is met, presenting a great challenge for the theory. 
However, on the basis of controlled experiments and com- 
prehensive analysis of the first principle model, we have 
demonstrated that the viscosity defines the intricate bal- 
ance between magnetic forces and hydrodynamic forces 
arising from the inertia of the particles and suspending 
liquid. The magnitude of these forces is inversely pro- 
portional to the viscosity, which can be independently 
controlled both in our experiment and theoretical model. 
We show that at a given frequency and amplitude of the 
energizing ac magnetic field, the viscosity of a suspending 
liquid controls the transition between snakes and asters; 
snakes emerge for smaller viscosities, while asters are fa- 
vored in more viscous liquids. 

In our experiments with liquid-air systems, the struc- 
tural transition from snakes to asters is consistently ob- 
served when the viscosity of suspending liquid is grad- 



ually increased. Our theoretical model, the fully non- 
linear Navier-Stokes equations coupled to the dynamics 
of individual magnetic particles, is reduced to a set of 
closed, time-averaged ordinary differential equations for 
particle positions and orientations interacting via mag- 
netic forces and effective mean hydrodynamic forces aris- 
ing due to oscillation of massive particles in a viscous 
liquid. The cause of these mean forces are Stokes drift 
and Rayleigh streaming. In contrast to the earlier study 
based on direct simulation of the Navier-Stokes equations 
[25|, here we obtain their analytic solutions, which pro- 
vide deep insight into the role of hydrodynamic flows, 
their detailed structure, and an overall balance of forces 
governing the self-assembly. The model is in good quali- 
tative (and some times quantitative) agreement with the 
experiments. 

Our experimental apparatus was similar to that de- 
scribed in Ref. [19[ . A ferromagnetic colloidal suspension 
was comprised of nickel microspheres with an average size 
of 90 /xm (Alfa Aesar Company). Due to defects in par- 
ticles, their magnetic moments are often strongly pinned 
and the particles behave as magnetically "hard" micro- 
spheres. The particles were dispersed at the liquid-air in- 
terface, where they were supported by a surface tension. 
To exclude the difference between the deep and shallow 
liquid layers, a circular glass beaker (5 cm in diameter) 
was filled with liquid depths of 5 cm and 5 mm. 

To vary the viscosity of the liquid, a range of water- 
sucrose solutions was prepared jlSj . The colloidal sus- 
pension was energized by an ac magnetic field, Hac = 
Hq sm{2TT ft) , with the frequency / and amplitude Hq = 
200 Oe, applied perpendicular to the interface. 

Selected experimental results arc summarized in Fig.[T] 
We observed the formation of magnetic snakes for values 
of the dynamic viscosity of the suspending liquid, 77, close 
to the viscosity of water. 77 w 1 mPas. With a gradual in- 
crease in 7], the snakes give way to asters, as illustrated in 
Fig. [TJ top panel. The transition is not sharp, it is asso- 
ciated with a wide transition region, as indicated by the 
error bars. Remarkably, the transition line is almost par- 
allel to the 77 axis above the viscosity of 77 « 1 — 2 mPa s. 
The bottom panel of Fig. [1] shows a characteristic time, 
Tf , for the formation of snake or aster as a function of i] 
for / = 40 Hz. After this time, the size of the developed 
structure almost did not change, the change of its rela- 
tive size was within 10%. Despite relatively large error 
bars. Tf gradually increases with the growth of rj. 

To obtain insights into the snake-aster transition, we 
significantly extend our model developed in Ref. 25| . We 
start with the model based on the fully nonlinear Navier- 
Stokes equation in the shallow water approximation, 

dth + V ■ (/iv) = , (1) 

c'tv+(v-V)v = z/ (V^v - av) - V/i + 7VV^/i 

+HosmiuJt)J2fir-r^)P^, (2) 
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FIG. 1: Top panel: Snake-aster phase diagram as a function 
of frequency, /, and viscosity, 77. The amplitude of the ac 
magnetic field is Hq — 200 Oe. Insets: Representative images 
of a snake (top) and an aster (bottom). Bottom panel: Char- 
acteristic time, Tf, for the formation of snakes (or asters) as 
a function 77 for / = 40 Hz and Hq = 200 Oe. For two points 
the vertical error bars are smaller than the markers. 



where v is the two-dimensional (2D), in-plane fluid ve- 
locity, h is the surface elevation, v is the kinematic vis- 
cosity, a is the friction with the bottom of the container, 
and 7 is the surface tension. The last term in Eq. ([2]) 
is representative of forces applied to the surface of the 
fluid through the particles, where Hq is the amplitude 
of the ac magnetic field, w is the frequency, the local- 
ized function / defines the shape of the particle, and 
Pi = {cos (f>i, sin (pi) is the orientation of the dipole mo- 
ment of the zth particle. In our study, we neglect the 
surface tension [26| and assume that / is given by delta 
functions (5(r— r^), which does not affect the basic physics 
of self-assembly but, more importantly, makes our model 
analytically tractable. The variables are scaled as fol- 
lows: coordinates r — > r/Zig, time t -> t^Jh^/g, veloc- 
ity V — >■ v/y/gTiQ, viscosity v — > v/hQ^/gho, where g is 
gravitational acceleration. In this dimcnsionlcss, rescalcd 
equation, a ^ ho ~ I. 

The motion of the particles on the surface of the fluid 
is described using Newton's equations 

mri + ixtYi = F,j + Titv - /3V/i , (3) 

I(i)i + fir<Pt = T, + kHo sin(wt)V/i X P, , (4) 

where m, /, ^t, IJr are the particle mass, moment of 
inertia, translational and rotational friction coefficients, 



Z^j^i ^ij 



respectively; /? = mg, Fj = J2j^i Fy and T^ 
are, respectively, the forces and torques on particle i due 
to magnetic and steric interactions with all other parti- 
cles, HtV is the Stokes' drag and —pVh is the movement 
along the surface gradient from gravity. The last term in 
Eq. ^ is the torque applied to each dipole moment in 
the direction of the projection of the vertical ac field on 
deformed surface [25|. 

In previous work, Eqs. ((l])-(|4]) were solved numer- 
ically to model snakes [25|. Here, we first analyti- 
cally find solutions of Eqs. ^ and ([2]) in an asymp- 
totic limit where we expand the surface deformation and 
liquid velocity with respect to the small parameter e, 
h = ho + ehi + e'^h2 + (e^) and v = evi -h e^vz + O {e') . 
The parameter e can be interpreted as the relative devi- 
ation of the locus /i(r, t) of the liquid-air interface from 
the equilibrium value h^. Moreover, by using the dimen- 
sionless viscosity v as a. small parameter, Eqs. ((T|) and 
(m were analytically solved up through the first order 
for the corresponding surface deformation and 2D veloc- 
ity fields induced by each particle individually to yield 
hi{r,t) = hr{v)e''^* + c.c. and Vi(r,t) = Vr(r)e'"* -fee, 
where c.c. denotes the complex conjugate. At the second 
order, time-averaged solutions h^ and V2 were sought and 
a corresponding analytic expression for V2, which detcr- 
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mines the mean flow, was obtained 

Using the explicit solutions /ii, /12, Vi, and V2 of the 
nonlinear Navier-Stokes equations, Eqs. ^ and ([2]), wc 
perform the time-averaging of Eqs. ([3]) and (j4]). As a 
result, we arrive at a closed system of ordinary differential 
equations for the particles in which all of the details of the 
complex hydrodynamic flows are effectively encapsulated 
in pairwise interactions 



mr, + ^jLtYi = ^ I F.y + Sj + fitv^2^^ - l3\'h 



j^i 






T,,-^V(/^W-/^(^■))xP. 



(5) 



(6) 



Here, the overline denotes complex conjugate and Sj = 

-2m[/3V|V/if^'^|2 + ^4{(v^^'.V)4''+c.c.}]/(a2 + m2w2) 

is the Stokes' drift term. To obtain the Stokes drift of 
each particle, we treated each term on the right-hand side 
of Eq. ^ independently. The last term in Eq. ((S]) is of 
much smaller order and can be neglected. 

Thus, in contrast to the earlier model [25[, where the 
dynamics of the particles is determined by Eqs. ([3]) and 
dl]) coupled to nonlinear equations ([T]) and ([2]), we sug- 
gest a much simpler and more transparent model in which 
the particle positions and orientations are described by 
Eqs. ([5]) and ([6]). Based on this model, we performed 
simulations with different numbers of particles ranging 
from 225 to 1000 [30[, with an initial configuration on a 






FIG. 2: Mean surface flows (V2) and first-order surface defor- 
mation {hr) induced by a single particle. Panels (a) and (b): 
Quadrupole streamlines produced by the potential (a) and ro- 
tational (b) components of the mean fiow, V2 . The magnitude 
of the velocity V2 decays as r~^. Panel (c): Composite of the 
potential and rotational flows. Panel (d): First-order surface 
deformation hr, which decays exponentially with r. The color 
in the image indicates the height. 



perturbed square lattice with a uniformly random orien- 
tation of the dipole moment and run on a GPU cluster. In 
addition to a significant reduction of computation time, 
roughly an order of magnitude speed up, the great advan- 
tage of our approach is gaining insight into the surface 
flows as the central ingredient underlying self-assembly. 

The overall analytic behavior of the mean surface flows, 
shown in Fig. [2j is similar to the large-scale quadrupolar 
flow seen from experiment. These flows are analo gou s 
to the mean flow produced by Rayleigh streaming |29| . 
The first-order flows (vi) arc time dependent, dipolar 
flows that oscillate in space and decay out exponentially, 
Vr(r) ex exp(— iA:r)/Y^, with fc w a; — ivki, ki = (w^ -|- 
a)/2; the behavior of hr{r) with r is similar to that of 
v,.(r), see Fig. [2jd). The second-order mean flow, V2, is 
time independent and is decomposed into the potential 
and rotational components, as shown in Figs. [2ja) and 
[^b), respectively. Both these counterparts have a long- 
ranged quadrupolar structure with a monotonic power- 
law decay oc r~^. The full mean flow is seen in Fig. [2Jc). 

Note that since the localized shape function was mod- 
eled by the delta function, the velocities and surface de- 
formations at all orders diverge at the center of each par- 
ticle. These divergences, however, have no effect on the 
system because each particle cannot influence itself and 
for each pair of particles, a short-range steric repulsion 



prevents them from getting close enough to feel the di- 
vergence. Earlier experimental data [20[ showed that a 
mean quadrupole flow was an essential ingredient for the 
assembly of snakes. Our model elucidates why this find- 
ing is true for both snakes and asters: Neither of these 
structures can be reproduced in simulations via Eqs. ([5]) 
and ([6]) unless the Stokes drift and the mean flow (i.e., 
the fields hr, v^, and V2) are properly determined. 

A critical test of the model is to recover the crossover 
in the behavior from snakes to asters that was seen ex- 
perimentally as a function of the liquid viscosity, 77, and 
the field frequency, /. The model successfully does so for 
a range of values of rj and /, see Fig. [3l top panel, where 
a = ho = p = 1. In qualitative agreement with the ex- 
periment, we observed snakes and asters formed for lower 
and higher values of i], respectively. Moreover, the depen- 
dence of time Tf for the formation of a dynamic structure 
on 77 exhibits a trend similar to the experimental one, see 
Fig. m bottom panel. Note that, in order to avoid depth 
dependence, the axes in Fig. [3] remain in dimensionless 
quantities. In the case where the viscosity was low and 
the frequency was high, the simulations yielded a clump- 
ing of particles primarily due to the lack of friction in the 
system. Alternately, when the viscosity was high and the 
frequency was low, the particles remained scattered due 
to overdamping and a lack of alignment along the changes 
in the surface height gradient. 

Because the mean flow induced by each particle has 
a long-range nature, it affects the dynamics of all other 
particles, leading to a highly nontrivial self-organization 
of the system. Figure |4] illustrates the formation of snakes 
and asters from an initially disordered distribution of par- 
ticles. The particles were dispersed uniformly inside a 
rectangle (snakes) or square (asters) with their magnetic 
moments oriented randomly. As Fig. |4] shows, asters and 
snakes are formed after a short transient (Tf), their orga- 
nization, e.g., anti-ferromagnetic order, closely resembles 
the experimental one. Starting from different initial con- 
ditions, e.g., square for the case of snake, often resulted 
in the formation of more than one snake or aster. 

In conclusion, we have demonstrated that the viscosity 
of the suspending liquid strongly affects the outcome of 
dynamic self-assembly and controls the structural tran- 
sition between self-assembled structures. Linear snakes 
are favored for small viscosities and circular asters for 
higher viscosities. Our novel model provides a nontriv- 
ial insight into how the large-scale mean flow - a non- 
linear effect caused by the strong coupling of oscillating 
particles with the initially equilibrium liquid - becomes 
a critical player that not only determines the shape and 
organization of the emergent dynamic structures but also 
keeps them stable. As we show, both inertia of the liq- 
uid and inertia of the particles are at the heart of the 
assembly process. 

While we were able to reproduce the main observed 
phenomenology, our method also has limitations. One 
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FIG. 3: Top panel: Snake-aster phase diagram as a function 
of the rescaled frequency and viscosity. Circles, triangles, and 
stars are, respectively for asters, snakes, and a mixture of seg- 
ments neither forming asters nor snakes. Inset: Comparison 
of experimental data (circles) and prediction of the model 
(solid lines) for ho = 1.5 mm and ho = 2 mm. Bottom panel: 
Formation time, Tf, of a structure as a function of fluid vis- 
cosity for f{g/ho)~^^'^ = 0.3. 



shortfall is that, being a completely two-dimensional 
model, the shallow water equations do not capture the 
liquid jets into the bulk produced by asters [2^ . However, 
our model is capable of replicating the overall structure of 
both snakes and asters. As snakes produce a largely 2D 
flow, it was natural to expect that they can be recovered 
from this model. Asters, however, produce a 3D toroidal 
flow and their appearance in this model was not expected. 
They form due to the propensity of the dipoles to align 
along the surface gradient for large characteristic wave- 
lengths, whereas entrainmcnt by the large-scale toroidal 
flow is important but less critical. Another reason for 
the lack of a good quantitative agreement between the 
experiment and model (such as in the inset of Fig. [31 top 
panel) is due to the neglected surface tension. Along with 
gravity, it presents another mechanism for surface wave 
generation but makes the model analytically intractable. 
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FIG. 4: Top row: Snake formation. Panel (a): Particles 
on a rectangular lattice with random orientation. Panel (b): 
Magnetic moments align along the surface gradient. Colors 
represent surface elevation h (red and blue show maxima and 
minima, respectively), arrows indicate particle magnetic mo- 
ments. Panel (c): Ferromagnetic chains are formed and are 
anti-ferromagnetically aligned, creating a snake. Bottom, row: 
Formation of an aster. Panel (d): Particles on a square lat- 
tice with random orientations. Panel (e): Magnetic moments 
align along the surface gradient. Panel (f): Ferromagnetically 
ordered chains are formed and an aster is assembled. 
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